Spin dynamics in HTSC cuprates: The singlet— correlated band (or t-J-V) model and 

its applications 
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So far calculations of the spin susceptibility in the superconducting state of cuprates have been 
performed in the framework of weak-coupling approximations. However, it is known that cuprates 
belong to Mott-Hubbard doped materials where electron correlations are important. In this paper 
an analytical expression for the spin susceptibility in the superconducting state of cuprates is de- 
rived within the singlet-correlated band model, which takes into account strong correlations. The 
expression of the spin susceptibility is evaluated using values for the hopping parameters adapted to 
measurements of the Fermi surface of the materials YBa2Cu307 and Bi2Sr2CaCu20s. We show that 
the available experimental data which are directly related to the spin susceptibility can be explained 
consistently within one set of model parameters for each material. These experiments include the 
magnetic resonance peak observed by inelastic neutron scattering and the temperature dependence 
of the NMR spin shift, spin-spin and spin-lattice relaxation rates in the superconducting state. 

PACS numbers: 74.25.Ha, 74.72.-h, 74.20.Mn, 74.25.Nf 



I. INTRODUCTION 

Derivations of theoretical expressions for the dynamical 
spin susceptibility of layered cuprates have been the fo- 
cus of many ongoing investigations, since various exper- 
imental quantities are directly related to the spin sus- 
ceptibility. Large amount of data sets exist on the tem- 
perature dependence of the Knight shift, spin-spin and 
spin-lattice relaxation rates as well as inelastic neutron 
scattering (INS) measurements. The most complete set 
of experimental data has been obtained for the YBaCuO 
compounds. 

The theoretical approaches to the spin susceptibility can 
be divided into two categories, the weak- and strong- 
coupling models. The former deals with a single-band 
Hubbard model with the effective Coulomb interaction 
U e taken to be of the order of the bandwidth in cuprates. 
Based on this assumption, the dynamical spin suscepti- 
bility can be calculated within a standard random phase 
approximation (RPA) approach. Extensive studies of 
NMR data were carried out within the framework of this 
model by Bulut and Scalapino et alm£££& and Mack et 
al£. In the meantime, the features observed by inelastic 
neutron scattering were studied theoretically by several 
groups^. However, in all of these calculations the ef- 
fect of strong electronic correlations were not taken into 
account properly. 

At the same time, there were a number of studies de- 
voted to analyzing the spin susceptibility within the 
strong-coupling t-J model, for which standard many- 
body perturbative methods do not work. For example 
the dynamical spin susceptibility was analyzed within 
the slave-boson approximation^, Mori-Zwanzig memory 
function formalism 10 - and with the Hubbard X-operators 



techniqu e) 11 ' 12 . It has been found that to a large extent, 
both weak and strong-coupling calculations give very 
similar results for the spin susceptibility. The respec- 
tive parameter values, however, differ drastically. Until 
now, there is no complete understanding whether both 
INS and NMR data can be explained consistently within 
one model and using the same parameter values of the 
given theory. 

In the present work we analyze this question in detail. 
Starting from the singlet-correlated band model, we use 
a well established decoupling procedure of the equations 
of motion and approximate higher order correlation func- 
tions so as to obtain an analytical expression for the 
dynamical spin susceptibility which takes into account 
strong correlations. It has been found by Hubbard and 
Jaini 3 " that the expression for the dynamical susceptibil- 
ity has a non-trivial correction due to strong correlation 
effects. Later Zavidonov and Brinkmannii have incor- 
porated an additional functional correction for the lower 
Hubbard sub-band (LHB) model, which accounts for lo- 
cal spin fluctuation effects. Both of these corrections are 
different from the conventional Pauli-Lindhard form and 
therefore they cannot be exactly included in the RPA 
approach. 

In the present paper we extend the previous analysis 
and present an analytical expression for the dynami- 
cal spin susceptibility in the upper Hubbard sub-band 
(UHB) in the superconducting state. We perform an ex- 
tended numerical evaluation of this analytical expression 
and find that most of the available experimental data 
which are directly related to the spin susceptibility can 
be explained consistently within one set of model param- 
eters. These experiments include the magnetic resonance 
peak observed by INS and the temperature dependence 
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of the NMR spin shift, spin-spin and the spin-lattice re- 
laxation rates, measured in the superconducting state. 
Note that in our analysis we restrict ourselves to op- 
timally doped high-temperature superconductors, since 
the pseudogap phenomenon cannot be explained within 
our model. Furthermore, in our analysis we take advan- 
tage of other available experiments, like the Fermi sur- 
face topology in various cuprate superconductors that 
is determined by high-resolution angle-resolved photoe- 
mission. Assuming a d. r 2_ J( 2-wave pairing symmetry, we 
propose an optimal set of parameters for the YBa2Cu30y 
and Bi2Sr2CaCu2C>8 compounds. 

This paper is organized as follows. In Sec.^we introduce 
the model system and present the analytical expression 
for the spin susceptibility in the superconducting state of 
cuprates. In Sec. IIIII and Sec.|W]we study the spin sus- 
ceptibility in the singlet-correlated band model by an- 
alyzing experiments in the superconducting state of the 
materials YBa2Cu3C>7 and Bi2Sr2CaCu20s. Summary 
and conclusions are given in Sec. Ivl 



II. DYNAMICAL SPIN SUSCEPTIBILITY IN 
THE SINGLET-CORRELATED BAND MODEL 



The starting point for our calculation is the model 
Hamiltonian 1214 
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where ipf d '' 7 (i^^' p j are composite copper-oxygen cre- 
ation (annihilation) operators of the copper-oxygen 
singlets states in the Cu02-plane. Furthermore Jy is 
the superexchange parameter of the copper spins (this 
coupling originates from the virtual hopping from LHB 
to UHB via the oxygen state). The number of doped 
holes is described by Si = ^f d,pd and Vy is an effective 
density-density interaction parameter. This parameter 
allows to account for the screened Coulomb repulsion 
and phonon (or plasmon) mediated interactions as well. 
Consequently, the last term is a Coulomb-like interaction 
between doped holes, which can be neglected, because it 
does not contribute to the spin susceptibility. 
The susceptibility is calculated from the general expres- 



X +- (q )W ) = -27ri((S+|Sr q ; 



(2) 



where the spin density operator S£ for the singlet- 
correlated band is written as 



(3) 



k+q- 



Here we neglected all the quasiparticle creation (annihi- 
lation) operators ?/>°' CT C&f'j corresponding to the LHB, 
because this band is assumed to be completely filled. 

The expression for the susceptibility is derived by the 
following procedure. First, we write down a complete 
set of equations of motion using the composite copper- 
oxygen creation (annihilation) operators tpf '' 

the copper-oxygen singlet states in the plane. Then, by 
means of a linear transformation we rearrange these equa- 
tions via Bogoliubov's quasiparticle operators into new 
sets of equations, which finally will be solved. An ex- 
pression for the susceptibility was previously derivecU* by 
utilizing the method of Heisenberg equations of motion 
in a small magnetic field. The advantage of the Green's 
function method is that it allows to obtain a formula for 
the susceptibility which contains both the itinerant (or 
quasi Fermi-liquid) part and the local spin fluctuation 
part in one general expression. 

The equation of motion for the relevant Green's function 
in the normal state (T>T C ) has been derived before by 
some of usi^. It is given by 
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where the factor P = (1 + 5)/ 2 is a doping dependent 
constant which arises due to the narrowing of the band 
in the so-called Hubbard-I approximation. 

In addition to Eq. (J2J it has been showni^ that 
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E^'-^ + q)((^Vai^: 



(5) 



Therefore, if we combine Eq. |0J and Eq. (JSJ) we get 

-^ktq|^q))= (6) 
-(^-^((-^^kfql^q 
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The equation of motion allows to derive the expres- 
sion of the dynamical spin susceptibility in the normal 
stated. 

For the superconducting state we need to perform Bo- 
goliubov's transformation 



(7) 



Consequently, the spin operator for the superconducting 
state will be written as 



(8) 
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Therefore in the superconducting state we need to 
construct additional equations for the Green's func- 

pd,l \,pd 
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Each of them has to be expressed via the -0^ ' operators, 
for example 
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and similar expressions can be 
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and 
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Equation i|l(J|) has the same 

form as in the normal state (Eq. ©), except that it is 
now adapted to be applied for the superconducting state. 
We note that in the conventional Fermi-liquid theory 
the anticommutator rule is given as c ka c k(T + c ka c k<? — 1- 
In the strong-coupling limit, however, this rule is 
modified^ due to the Coulomb repulsion. For this rea- 
son we have abbreviated the terms which are present in 
the conventional weak-coupling Fermi-liquid approach 
in the superconducting state by the truncated Green's 



function 



V"k ^k+q' 

terms on the right hand side of Eq. are due to the 
spin modulation S£- 

With the help of these equations we are able to construct 
the equations of motion which are needed to calculate the 
spin susceptibility in the superconducting state. The first 
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one is given as 



(uj-E p + E k )((oi d ' l al pd \SZ 



(11) 



2vr 



(M k Mp + v k v p ) (rip - n k ) 
(u k u p + v k v p ) {( J q - t p ) 



N 

-(J q -ik)n k }((5+|5Z q )) 

- (Pu k u p + (P- l)v k v p ) jL «5+ |5Z q » , 



"k+q 



s: 



y k+q |"-q 



5': 



-w k u k+q (Ul^' 1 ^ 

+^k + q((^if^|^ q )) 



+ ^k + q((V^V_tiq|^q 



Doing so we get 
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and similar expressions occur for / / a]^e£?^ \S_ q 
(alt^lSZ^)) and ((<'V_^ |SZ q )). Fur- 

are the 



thermore, n k = u k 

occupation numbers in the superconducting state. We 
further make use of the identity 
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With help of this relation the susceptibility is calculated 

as 
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where the functions x k 



(13) 



xt' (q,w) 



1 + JqXo (q,w)+n(q,w) + Z(q,w) 



where the superexchange interaction between the cop- 
per spins is J q = J\ (cosq x + cosq y ), with J\ being the 
superexchange interaction parameter between the near- 
est neighbour copper spins. The function Xo~'~ "(q> a; ) 
is a BCS-like susceptibility and II(q, oS) is a function 
which results from strong correlation effects and has been 
determined^ before. It is given by 



P 
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The function Z(q, w) has its origin in the fast fluctuation 
of the localized spins and it is calculated as 
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w k = l(l-£k/Ek) and z k = u k v k = A k /(2E k ) are 
the conventional coherence factors. Furthermore T is 
an artificially introduced damping constant and E k = 



(e k — //) + A£ is the energy of Bogoliubov's quasipar- 
ticles in the superconducting state. The energy disper- 
sion in tight-binding approximation for a quadratic two- 
dimensional lattice is given as 



P [2h (c 



(16) 
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(cos k x + cos k y ) 
+ 4< 2 (cos k x cos k y ) 
+ 2i 3 (cos 2k x + cos 2k y ) 
+ 2<4 (cos 2k x cos k y + cos 2k y cos k x ) 
+ 4i 5 (cos 2k x cos 2ky)\ , 



where the model parameters ti,t2— correspond to 
nearest-neighbour (NN), next-nearest-neighbour 
(NNN), and further distant hopping, respectively. For 
simplicity we do not consider hopping between layers. 
Further we note that at optimal doping the number of 
doped holes per unit cell in one Cu02-layer is 0.165. In 
bilayer compounds therefore we have S = 0.33 with a 
corresponding factor P = (1 + S)/2 ~ 0.7 near optimal 
doping. 

The mechanism that causes the pairing in cuprates is 
still being debated and the origin of the interactions de- 
scribed by Vij in Eq. Q arc unknown. Therefore we 
have introduced the superconducting gap function Ak 
phenomenologically into our model. Assuming a d x 2_ y 2 
pairing symmetry it is given by 



A k (T) = ^ (cos k x - cos ky) tanh (l.76y/T c /T - 1 

(17) 

where Ao is considered to be a model parameter. We 
would like to point out that this formula is a fit to the 
solution of the Eliashberg strong-coupling gap equation. 

In the forthcoming sections we will analyze sev- 
eral experiments in the superconducting state of 
the materials YBa2Cu307 and Bi 2 Sr 2 CaCu208. Wc 
would like to summarize at this point the pa- 
rameters which we used to perform our analysis. 
The tight-binding hopping parameters {Pt\...Pt^) are 
adopted from fits to the measured Fermi surfaces 
of these two materials. For YBa2Cu30y we hav©i£ 
(fi,Ph...Pt 5 ) = (119,147,-36.5,-2.4,32.4,-1.8) meV, 
while for Bi 2 Sr 2 CaCu20 8 we takei£ (n,Pt 1 ...Pt 5 ) = 
(49.4, 73.9, -12.0, 16.3, 6.3, -11.7) meV. The correspond- 
ing Fermi surfaces are shown in Fig. ^ Note that they 
are quite different for these two materials. 

Other model parameters include the gap parameter 
which is assumed to be in the order of Ao ~ 10 — 30 meV 
and the superexchange interaction parameter of the cop- 
per spins which is J\ ~ 100 — 140 meV. In our analysis 
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Figure 1: Fermi surfaces of YBa2Cu3C>7 (dashed) and 
Bi2Sr2CaCu2C>8 (solid) according to photoemission experi- 
ments. 

we proceed as follows. We assume the Fermi surface as 
given from fits to photoemission data and thus fix the val- 
ues of the hopping parameters for both materials. Based 
on this assumption we analyze neutron scattering exper- 
iments which allows us to determine the values of the 
model parameters Ao and J\. Then we move to NMR 
experiments and calculate the temperature dependences 
of spin shift, spin-spin relaxation and spin-lattice relax- 
ation rates, utilizing the parameter values obtained be- 
fore. Thus in our analysis we try to establish a connection 
between the different experiments. 



III. NEUTRON SCATTERING ANALYSIS 




Figure 2: Calculated frequency and momentum dependence 
of Imx + '~ (q, for YBa2Cu307 (top) and Bi2Sr2CaCu20s 
(bottom). 



Magnetic inelastic neutron scattering experiments probe 
directly the imaginary part of the dynamical spin suscep- 
tibility Imx ~ (q, w). The experiments indicate a sharp 
resonance in the magnetic excitation spectrum of opti- 
mally doped YBa 2 Cu 3 OTi2i2a and B^SraCaCuaOgSiiS 2 . 
compounds at a frequency lu ~ 41 meV, near the antifer- 
romagnetic wave vector Q = (tt,tt). Consequently, there 
should be a large peak in the imaginary part of the spin 
susceptibility Imx + '~ (Q,t^) at the same frequency. In 
the conventional weak-coupling scenario this feature was 
studied extensively by various authors^, who connected 
the appearance of the resonance peak to a collective 
spin-density wave mode formation. Special experimen- 
tal features like the effect of orthorhombic distortions 23 
and bilayer splitting 2 ^ on the magnetic excitations were 
also studied theoretically 2 ^* 2 ^ within the weak-coupling 
model. In the strong-coupling limit previous calculations 
were carried out by some of us^ and the dependence 
of the position of the resonance peak on the model pa- 
rameters was studied extensively. We will not repeat 
these considerations here. Note only that within our 



model the position of the neutron scattering resonance 
peak is determined mainly by the magnitude of the su- 
perconducting gap Ao and the superexchange parame- 
ter J\. In particular, the superconducting gap parame- 
ter Ao determines the size of the transparency window 
in Imx +, ~ (q, w) which is approximately u> ~ 2Ao. In 
this region a sharp delta-like peak appears in the imag- 
inary part of the susceptibility if the resonance condi- 
tion 1 + JqReXo '~{<i,uj) + ReU(q,uj) + ReZ(q,u) = 
(see Eq. (|13|l ) is fulfilled. Our calculations indicate that 
for YBa 2 Cu307 (Bi 2 Sr 2 CaCu 2 0s) the value of the gap 
parameter should be Ao = 24 meV (25 meV). The cor- 
responding values of the superexchange parameter are 
determined as J\ = 90 meV (110 meV). For these val- 
ues the resonance condition is fulfilled and a clear peak 
appears in the imaginary part of the susceptibility near 
uj ~ 41 meV, for both materials. In Fig. [3 we display the 
calculated momentum and frequency dependence of the 
imaginary part of the susceptibility. Note that the height 
of the resonance peak depends on the artificially intro- 
duced quasiparticle damping T, therefore the values of 
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model calculations match the experimental observation 28 
of incommensurability. For the Bi 2 Sr 2 CaCu20g com- 
pound, however, our results indicate that the incom- 
mensurability of the magnetic excitations is much weaker 
than in YBa2Cu3C>7 due to the difference in the Fermi 
surface topology. This result could be tested by further 
experiments. 



IV. NMR ANALYSIS 



A. Knight shift 

.5 

In order to calculate the Knight shift in the supercon- 
| ducting state we need to calculate the susceptibility in 

the limit q — > 0, u> = 0. The BCS susceptibility 
Xo~ (Q; lu ) converts to the Yosida result 



Figure 3: Intensity plot of the imaginary part of the sus- 
ceptibility Imx + '~ (q, w) near Q for YBa2Cu307 (top) and 
Bi2Sr2CaCu208 (bottom) for uj — 30 meV. 



The functions II(q, uS) and Z{(\ 1 lu) are approximated as 



il(q-»0,w = 0) 



^/(S k ) (19) 
df (E k ) 



N 



k 

P(3 

k 

P ~ ~N~ 



df (E k ) 



Imx + '~ (q, u>) in Fig. [21 are arbitrary. Furthermore, the 
experimentally reported^ downward dispersion branch 
for YBa2Cu3C>7 with respect to lu is reproduced by our 
model calculations as can be seen from an inspection 
of Fig. |21 (top). We would like to point out that for 
the Bi2Sr2CaCu208 compound we do not find a similar 
dispersion branch. There are, however, no experiments 
available which would allow a comparison. 

Let us now turn to the examination of magnetic excita- 
tions at lower frequency, where measurements^ 8 , in the 
YBaCuO compounds indicate well defined incommensu- 
rability in the magnetic excitation spectrum. For the 
material Bi2Sr2CaCu20s, however, the data setsSS are 
inconclusive due to experimental difficulty. Here we re- 
port significant differences in the low frequency excita- 
tions for the materials YBa2Cu3C>7 and Bi2Sr2CaCu20s. 
Similar conclusions were reached previously by Norman^ 
in the conventional weak-coupling scenario. In Fig. we 
show an intensity plot of the imaginary part of the sus- 
ceptibility around the antiferromagnetic wave vector Q 
calculated for lu = 30 meV, for both materials. By exam- 
ination of the figure we conclude that for YBa2Cu3C>7 the 



and 



Z(q -> 0,lu = 0) ~ P. 



(20) 



In the long-wave limit therefore the susceptibility is given 
by the simple expression 



x +>-(q^0,w = 0) = 



Xp 



1+P + 5/P+ (2Ji - fi/P) X p ' 

(21) 

With the help of this relation the spin shift can be cal- 
culated according to 



K s cx A X + '-(q^Q,LU = Q), 



(22) 



where A represents the appropriate hyperfine coupling 
constant. Note that this expression refers to the spin con- 
tribution to the magnetic shift. In addition there is an 
orbital (chemical) shift which, however, is independent 
of the temperature. We will calculate the temperature 
dependence of the normalized spin shifts K S (T)/ K S (T C ). 
In this way the hyperfine coupling constants in Eq. I|22|) 
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Figure 4: Temperature dependence of the reduced spin shift Figure 5: Temperature dependence of the reduced spin shift 
in YBa2Cu307 for no exchange interaction (Ji = eV). The in YBa2Cu307. The energy gap is Ao = 24 meV. The exper- 
experimental points are taken from Barrett et all— . imental points are taken from Barrett et all— . 



cancel out, simplifying our analysis. Furthermore, note 
that the quasiparticle damping is Y — > + . In Fig. 0] 
we display the calculated temperature dependence of the 
spin shifts for YBa 2 Cu307, along with the experimental 
points of Barrett et al^, for vanishing super-exchange 
interaction J\ = eV and different values of the gap pa- 
rameter Ao. We observe that below T c the spin shifts 
depend strongly on the magnitude of the gap parame- 
ter Ao. This behavior has also been found for the RPA 
susceptibilit}^. 

Next we consider how the spin shift depends on the su- 
perexchange interaction J\. In Fig. [5] the calculated spin 
shifts for different values of the superexchange interac- 
tion parameter Ji are shown. We see that the tempera- 
ture dependence of the Knight shift does not significantly 
change by adjusting the parameter J\, Also, contrary to 
the RPA scenario, the superexchange coupling J\ reduces 
the rapid decrease of the Knight shift. By analysis of the 
figure we conclude that the optimal set of parameters to 
describe the experimentally observed temperature depen- 
dence of the spin shift for YBa2Cu3C>7 is Ao = 24 meV 
and Ji = 90 meV for the given Fermi surface. These val- 
ues are in perfect agreement with those determined from 
the fit to neutron scattering experiments in the previous 
section. 

Let us now turn to the examination of the spin shift in 
the Bi 2 Sr 2 CaCu208 compound. Experimental results in- 
dicate a similar behaviour as in the YBa2Cu3C>7 material: 
the spin shift decreases rapidly upon entering the super- 
conducting state. The calculated spin shifts also show a 
similar dependence on the model parameters Ao and J\ . 
We will not repeat the analysis of these dependences and 
show instead in Fig. the final result of our calculations 
for the spin shift in Bi2Sr2CaCu20s along with the exper- 
imental points of Ishida et al& The parameters used for 
the calculation are Ao = 24 meV and J\ = 110 meV. We 
notice that again these values almost coincide with those 
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Figure 6: Temperature dependence of the reduced spin shift 
in Bi2Sr2CaCu2C>8. The experimental points are taken from 
Ishida et aliS. 



determined by the analysis of neutron scattering experi- 
ments in this compound. By examination of the figure we 
see that the calculated temperature dependence of spin 
shift gives a satisfactory fit to the experimental data. 

Finally we conclude that it is possible to account for both 
the neutron scattering resonance peak and the tempera- 
ture dependence of the spin shift in the superconducting 
state of YBa2Cu3C>7 and Bi2Sr2CaCu20s consistently 
within the same set of parameters for each material. 
Next, we calculate the temperature dependence of dy- 
namical NMR quantities, the spin-spin and spin-lattice 
relaxation rates. 
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Figure 7: Temperature dependence of the spin-spin relaxation 
rate in YBa2Cu307 for no exchange interaction (Ji = eV). 
The experimental points are taken from Stern et alf— . 



Figure 8: Temperature dependence of the spin-spin relaxation 
rate in YBa2Cu307. The energy gap is Ao = 22 meV. The 
experimental points are taken from Stern et alf— . 



B. Spin— spin relaxation 

The nuclear spin-spin relaxation rate is calculated from 
the expression 3 ^ 



T 2G K 



^ 63 F,|(q) 2 ( J R eX + '-(q,c = 0)) 2 (23) 



where 63 F||(q) = [A\\ + 2B (cos(<f x ) + cos(q y ))] 2 is the 
hyperfine form factor. The values of the hyperfine cou- 
pling constants are taken as B ~ 0.4 /zeV and A\\ ~ — 4B. 
The spin-spin relaxation rate allows us to study the real 
part of the susceptibility i?ex +, ~(q, u = 0) near the an- 
tiferromagnetic wave vector Q. We will calculate the 
temperature dependence of the normalized spin-spin re- 
laxation rate T 2 ^ 1 (r)/T 2 ^ 1 (r c ) in the same way as we 
analyzed the Knight shifts. Note that for the evaluation 
of the spin-spin relaxation, the real part of the suscepti- 
bility is calculated by taking the quasiparticle damping 
r — > + . Otherwise, due to the behavior of the coherence 
factors, a large increase in the spin-spin relaxation rate 
occurs near T c upon entering the superconducting state 
as has been discussed in Refi 3 ^. 

In Fig. we display the calculated spin-spin relaxation 
rates for YBa2Cri307 along with experimental points 
from Stern et al^ for no superexchange interaction 
J\ = eV. We observe that the results show a simi- 
lar temperature dependence as in the RPA approach 6 . 
Generally, the temperature dependence of the spin-spin 
relaxation rate is less sensitive to the change of the gap 
parameter than the spin shift. We see that for the hy- 
pothetical case of no interaction we can account for the 
observed temperature dependence of the spin-spin relax- 



ation rate. Next we wish to study the behavior of the 
spin-spin relaxation for different values of the superex- 
change interaction parameter J\. 

In Fig. |S| the temperature dependence of the spin-spin 
relaxation rate is shown for various values of the superex- 
change interaction parameter J\. We see that we get a 
reasonable agreement with the data using the parameter 
values Ao = 22 meV and J\ = 90 meV. The magnitudes 
of these parameters are in agreement with those obtained 
by the analysis of neutron scattering and NMR spin shift 
experiments for YBa2Cu30y . 

Concerning the I^Si^CaC^Og compound our calcula- 
tions indicate a very similar behaviour as in YBa2Cu30y, 
if we utilize the parameter values from the previous sec- 
tions Ao = 24 meV and J\ — 110 meV. Unfortunately 
the spin-spin relaxation rate has not yet been measured 
in Bi2Sr2CaCu20g , thus we have no basis for comparison 
with experiments. 



C. Spin— lattice relaxation 

The nuclear spin-lattice relaxation rate is calculated ac- 
cording to the expression 3 ^ 



T ip oc — 



1,0' 



lim 



Imx +t (q,w) 



(24) 



where (3 denotes the field direction and j3' are the direc- 
tions orthogonal to the field. Furthermore a designates 
the nucleus under consideration. 

In order to calculate the imaginary part of the spin sus- 
ceptibility irox + '~(q, ijJ — > 0) we introduced a finite 
quasiparticle broadening Y — 3fcsT c ~ 2 meV, follow- 
ing the analysis of Bulut and Scalapinc 5 . Furthermore, 
the form factors in Eq. I|24|l are given by 
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Figure 9: Temperature dependence of the spin-lattice relax- 
ation rate in YBa2Cu3C>7 for no exchange interaction (Ji = 
eV). The experimental points are taken from Takigawa et 



63 i^(q) = [A p + 2B {cos(q x ) + cos(q y ))} 2 , (25) 
17 F^(q) = 2(Clcos 2 (q x /2) + Clco S 2 (q y /2)). 

The values of the hyperfine coupling constants are taken 
as B ~ 0.4 fjjeV, A {{ ~ -4B, A ± ~ 0.755, C\\ ~ 0.6B, 
and Cj. ~ 0.32B. 

In Fig. [5] we display the calculated spin-lattice relaxation 
rates for YBa2Cu3C>7 along with the experimental points 
of Takigawa et al>21, when the superexchange interaction 
is J\ = eV. We observe that the temperature depen- 
dence varies strongly when adjusting the gap parameter 
Ao- As for the spin shift and spin-spin relaxation rate 
calculations it is possible to fit the experimental data 
even without taking into account the interaction. 

Next we consider the effect of the superexchange parame- 
ter J\. In Fig.m|the spin-lattice relaxation rate is shown 
for different values of J\. We note that upon chang- 
ing the values of J\ the spin-lattice relaxation rate T^} 
changes the same way as it does in the RPA case if the 
parameter value of the effective Coulomb interaction U e 
is changed. Namely, the parameter J\ has no significant 
impact on the temperature dependence of the spin-lattice 
relaxation rate in the superconducting state. Upon fur- 
ther examination of the figure we see that we get a rea- 
sonable agreement with experimental observation using 
the parameter values Ao = 22 meV and J\ = 90 meV. 
These parameters agree with those we determined before 
in the previous Sections. 

Next we examine the spin-lattice relaxation rate in the 
Bi2Sr2CaCu208 compound. Upon changing the model 
parameters Ao and J\ the spin-lattice relaxation behaves 
much the same way as in YBa2Cu307. We show in Fig. 
ITTl the final result of our calculations of the spin-lattice 
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Figure 10: Temperature dependence of the spin-lattice relax- 
ation rate in YBa2Cu307. The energy gap is Ao = 22 meV. 
The experimental points are taken from Takigawa et alp—. 
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Figure 11: Temperature dependence of the spin-lattice re- 
laxation rate in Bi2Sr2CaCu2 0s. The experimental points 
are taken from Ishida et al,' 3 (squares) and Takigawa et all— 
(circles). 



relaxation rate in Bi2Sr2CaCu20s, along with the exper- 
imental points of Ishida et ali^ (squares) and Takigawa 
et aliSi (circles). The parameters used for the calcula- 
tion are Ao = 23 meV and J\ = 110 meV. Note that the 
values of these parameters are again almost the same as 
those that we used before when we analyzed neutron scat- 
tering and spin shift experiments. A particularly interest- 
ing feature can be found when comparing the spin-lattice 
relaxation rates in YBa2Cu307 and Bi2Sr2CaCu20s at 
low temperatures. A close inspection of the correspond- 
ing figures (Fig. EI] and Fig. fTT|) shows that in the for- 
mer case 63 T 1 ~ 1 (T) practically vanishes at temperatures 
T< 20 K, while in the latter case the relaxation rate 
seems to vanish only at T~ K. Note that both of these 
dependences are reproduced by the model calculations. 
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Figure 12: Temperature dependence of the anisotropy ra- 
tios 63 T 1 ~ 1 / 17 T 1 ~ 1 and ^T^/ 63 !^ 1 in YBa 2 Cu 3 7 (squares). 
Data for YBa2Cu4 08 are also shown for comparison (circles). 



We are also interested in the anisotropy ra- 
tios 63 T 1 ^/ 63 T 1 ; 1 and 63 / 17 T^ measured in 
YBa2Cu30y. Experimental evidenca^MSdi points 
toward a field dependence of these quantities. Our 
theoretical results have to be compared with data in 
zero or small external fields. We display the calculated 
anisotropy ratios in Fig. 1121 The experimental points 
for 63 T , 1 " C 1 / 17 T 1 ; 1 are from Martindale et ali22, whereas 
those for 63 T 1 -^/ 63 T 1 ; 1 are taken from Takigawa et 
al41 (squares). For comparison we also plotted the 
anisotropy in the YBa 2 Cu 4 08 compound (circles) 
measured by Bankay et al4i. For the calculation we 
used the same parameters as before, Ao = 22 meV 
and Ji = 90 meV. We see that we can account for the 
anisotropy ratio 63 7 1 1 ^J ) / 63 T 1 ~ 1 , but it is not possible 
to reproduce the temperature dependence of the ratio 
63 ^lc / 17 ^ie • Our results for the weak-field anisotropy 
ratio 63 r i "] ) / 63 T 1 - 1 for YBa 2 Cu 3 7 agree with the 
weak-coupling calculations of Bulut and Scalapinc-, 
which are based on a square Fermi surface with nearest 
neighbour hopping only. However, our calculations 
disagree with the analysis of Mack et al£, where the 



Fermi surface for YBa 2 Cu30y was assumed to be quite 
different (similar to that which we used here to describe 
the Bi 2 Sr 2 CaCu208 compound). Therefore we conclude 
that the Fermi surface topology plays an important 
role in the description of the weak-field anisotropy 
ratio 63 T i ;] ) / 63 r i ~ 1 . Note that the experimentally 
reported weak-field ratio 63 T , 1 ~ 1 / 17 T 1 ~ 1 could also not 
be reproduced within previous weak-coupling RPA 
calculations^. 



V. CONCLUSIONS 

In summary we have determined the spin susceptibility in 
cuprates within a special Hubbard model which includes 
strong correlation effects. It has been found that the sus- 
ceptibility in the strong-coupling limit is different from 
the standard Pauli-Lindhard formula. In particular, two 
correction functions were determined in the supercon- 
ducting state. The first one, n(q, us), found originally by 
Hubbard and Jaim— in the normal state, originates from 
the anticommutator rule which is modified due to the 
Coulomb repulsion, whereas the function Z(q, to) has its 
origin in the fast fluctuations of the localized spins and 
was previously discussed by Zavidonov and Brinkmann 11 
for the normal state. 

We analyzed inelastic neutron scattering and NMR 
data in the superconducting state of the opti- 
mally doped high-T c superconductors YBa 2 Cu307 and 
Bi 2 Sr 2 CaCu 2 08. In our analysis we have taken into ac- 
count the experimentally measured topology of the Fermi 
surface, which is quite different for these two materials. 
We found that on the whole the results within the strong- 
coupling and weak-coupling limits agree with each other. 
Based on the results of our numerical calculations we 
conclude that strong correlation effects, i.e. the effect of 
the functions n(q, u>), Z(q,uj) on the susceptibility can 
be modeled in the weak-coupling approach by an appro- 
priate redefinition^ of the effective Coulomb interaction 
parameter U e . In particular, the non-physical value of 
U e in the weak-coupling limit (sometimes U e > t) be- 
comes understandable. In terms of our model this can be 
explained by the impact of the n(q, uj) function, which 
is of the order of tx^ ~ (q, uj). As concerns the Z(q,uj) 
function, our calculations show that if it is not included, 
the numerical values of /?7ix + '~(q, uj) can for some set of 
parameters become negative, implying the importance of 
this correction. 

In the framework of the singlet-correlated band model 
we found it possible to describe the available experi- 
mental data in the optimally doped YBa 2 Cu3 07 and 
Bi 2 Sr 2 CaCu 2 08 compounds within one set of model pa- 
rameters for each material. These optimal sets of param- 
eters are given by A = 23 meV (±5%), J\ = 90 meV for 
YBa 2 Cu 3 7 and A = 24 meV (±5%), J x = 110 meV 
for Bi 2 Sr 2 CaCu 2 0s. The only experiment which could 
not be reproduced is the temperature dependence of the 



11 



weak-field anisotropy ratio 63 T 1 ; 1 / 17 ^ 1 in YBa 2 Cu 3 7 . 
As concerns the Bi 2 Sr 2 CaCu20g material, more experi- 
mental data would be desirable to test our model further, 
for example low energy inelastic neutron scattering and 
spin-spin relaxation rate measurements. 
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